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We study the prospects to measure the CP-sensitive triple-product asym- 
metries in neutralino production e + e~ — > Xi Xx an d subsequent leptonic two- 
O body decays Xi ~^ £r£> @R ~^ Xi^; f° r @ = e )M; within the Minimal Su- 

CN ■ persymmetric Standard Model. We include a full detector simulation of the 

. International Large Detector for the International Linear Collider. The simu- 

lation was performed at a center-of-mass energy of yfs = 500 GeV, including 
' the relevant Standard Model background processes, a realistic beam energy 

spectrum, beam backgrounds and a beam polarization of 80% and —60% for 
the electron and positron beams, respectively. In order to effectively disentan- 
gle different signal samples and reduce SM and SUSY backgrounds we apply 
a method of kinematic reconstruction. Assuming an integrated luminosity 
of 500 fb" 1 collected by the experiment and the performance of the current 
ILD detector, we arrive at a relative measurement accuracy of 10% for the 
CP-sensitive asymmetry in our scenario. We demonstrate that our method 
of signal selection using kinematic reconstruction can be applied to a broad 
class of scenarios and it allows disentangling processes with similar kinematic 
properties. 
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1 Introduction 



Supersymmetry (SUSY) [HI2] is one of the best motivated candidates for physics 
beyond the Standard Model (SM). Besides providing a unification of the strong and 
electroweak gauge couplings and a suitable cold dark matter candidate, the lightest 
stable SUSY particle, SUSY could offer new sources of CP violation [3H6]. In the 
Minimal Supersymmetric Standard Model (MSSM), the complex parameters are 
conventionally chosen to be the Higgsino mass parameter /x, the U(l) and SU(3) 
gaugino mass parameters M\ and M3, respectively, and the trilinear scalar coupling 
parameters Af of the third generation sfermions (/ = b,t,r), 

H = |/x|e^, Mi = iMile** 1 , M 3 = |M 3 |e^ 3 , A f = \A f \e^ A f . (1) 

The sizes of these phases are constrained by experimental bounds from the electric 
dipole moments (EDMs). Such experimental limits generally restrict the CP phases 
to be small, in particular the phase M [7j. However, the extent to which the EDMs 
can constrain the SUSY phases depends strongly on the considered model and its 
parameters [7HT3]. 

Due to cancellations among different contributions to the EDMs, large CP phases 
can give CP- violating signals at colliders, as shown for example in Ref. [13]. It is 
important to search for these signals, since the cancellations could be a consequence 
of an unknown underlying structure that correlates the phases. In addition, the 
existing EDM bounds could also be fulfilled by including lepton flavor violating 
couplings in the slepton sector [9]. 

Thus, direct measurements of SUSY CP-sensitive observables are necessary to 
determine or constrain the phases independently of EDM measurements. The phases 
can change SUSY particle masses, their cross sections, branching ratios [HHT8] , and 
longitudinal polarizations of final fermions [TH] . Although such CP-even observables 
can be very sensitive to the CP phases, CP-odd (T-odd) observables have to be 
measured for a direct evidence of CP violation. 

CP asymmetries in particle decay chains can be defined with triple products of 
final particle momenta [201121] . Due to spin correlations, such asymmetries show 
unique hints for CP phases already at tree level. Thus, triple product asymme- 
tries have been proposed in many theoretical papers. For the Large Hadron Col- 
lider (LHC), triple product asymmetries have been studied for the decays of neu- 
tralinos [22H23], stops [T31l2oW]. sbottoms [2gl[22|, and staus [3D]. In a Monte 
Carlo (MC) analysis for stop decays [2nj[27], it could be shown that the decay 
chain can be reconstructed and asymmetries be measured at a 3a level for luminosi- 
ties of the order of 300 fb _1 . At the International Linear Collider (LLC) [311431] a 
clearer identification and a more precise measurement is expected to be achievable. 
However, in this context only theoretically-based papers exist: for instance, neu- 
tralino production with two- [2"2"jl3"oTH2"] and three-body decays [I3THT] . charginos 
with two- [HH53] and three-body decays [1131 155] . also with transversely polarized 
beams [56H6T] . have been studied. 
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Figure 1: Schematic picture of neutralino production and decay. 

Therefore, we present in this paper the first experimentally-oriented analysis 
with regard to the observation of CP asymmetries on the basis of a full detector 
simulation. To show the feasibility of a measurement of triple product asymmetries, 
we focus on neutralino production [3"5| l3"9] 

e + + er -+ t + xl (2) 

with longitudinally polarized beams and the subsequent leptonic two-body decay of 
one of the neutralinos into the near lepton 

x°i ->• e R +t N , (3) 

and that of the slepton into a far lepton 

Ir Xi + ^S £ = e lt M. (4) 

Fig. [T] shows a schematic picture of neutralino production and its decay chain. The 
CP-sensitive spin correlations of the neutralino in its production process allow us to 
probe the phase of the Higgsino mass parameter /x, and the gaugino parameter M\. 

In order to effectively disentangle different signal samples and reduce SM and 
SUSY backgrounds we apply a method of kinematic reconstruction. A similar ap- 
proach has been studied successfully for the LHC [27] . The kinematic reconstruction 
was also considered to study selectron and neutralino properties at the ILC [52H53] . 
Here, we demonstrate that it can be used as an effective signal selection method, 
greatly improving the sensitivity to the effects of CP violation. In particular, com- 
pared to the previous studies of process (J2J, Ref. [39], we are able to suppress slepton 
and WW contamination to 0(10%) level. 

To investigate in detail the prospect of measuring CP-sensitive observables at 
the ILC we perform a full detector simulation of the International Large Detec- 
tor (ILD) concept. We include all relevant SM and SUSY background processes 
in our study, simulated with a realistic beam energy spectrum (beamstrahlung 
and initial state radiation (ISR)), beam backgrounds and a beam polarization of 
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(P e -,P e +) = (0.8, —0.6), which enhances the cross section of our signal and the size 
of the asymmetry. We apply the method of kinematic reconstruction to a preselected 
sample of signal event candidates in order to efficiently reject any background and 
to disentangle the decays of x® from x®- We determine the CP asymmetries with 
the selected signal events and study the sensitivity to determine the values of the 
CP phases via a fit. 

The paper is organized as follows. In Sec. H]we introduce the theoretical frame- 
work for the used CP-sensitive observable and we apply it for the studied process in 
the chosen benchmark scenario. Section [3] discusses the kinematic aspects of signal 
versus background selection. Section H] treats the full detector simulation. Finally 
in Sec. the SUSY parameters including the CP phases are derived via a fit of the 
CP-odd asymmetries together with masses and cross sections. Appendix R1 provides 
details for the reconstruction of W and I production and App. [B] recapitulates the 
neutralino mixing and its parameters. 

2 Theoretical framework 

2.1 CP-odd observables 

CP-violating observables in collider-based experiments are based on the invariance 
under CPTn, where C is charge conjugation, P stands for parity transformation and 
T for time reversal. The index N denotes 'naive' time reversal, i.e. time reversal 
but without interchanging initial with final states and therefore can be tested in 
collider-based experiments. At tree level of perturbation theory, observables odd 
under Tn transformations are also odd under the 'true' time reversal T. 

Therefore, it is useful to categorize CP- violating observables in two classes [65J: 
those that are even under Tn and those that are odd under Tn operation. Under 
the absence of final state interactions (FSI), CPTN-even operators relate TN-odd 
symmetries uniquely with CP-odd transformations [20]. Contrary, CPTN-odd op- 
erators (i.e. CP-odd but TN-even) can have nonzero expectation values only if FSI 
are present that give a non-trivial phase (absorptive phase) to the amplitude. Such 
a phase can arise for instance in loop diagrams. 

Examples for TN-odd observables are triple products that arise from the terms 
e [Pi;P2;P3;P4] ; where Pi are 4- vectors representing spins or momenta and e is the 
antisymmetric Levi-Civita tensor. Consequently, such TN-odd signals can only be 
observed in processes where at least four independent momenta (or their spin ori- 
entations) are involved. The e-tensor can then be expanded in a series of four triple 
products e[pi,P2,P3,P4] = E1P2 ' (P3 x P4) ± • • • that can be evaluated in a suitable 
and specific kinematical system. The TN-odd asymmetries are then composed by 
the corresponding triple products. 

2.2 Neutralino production and decay processes 

Neutralinos are mixed states of the supersymmetric partners of the neutral gauge 
and Higgs bosons and depend on the phases <f>\ and see App. [B] CP- violating 
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effects in the neutralino production and decay arise at tree level and can lead to 
CP-sensitive asymmetries due to neutralino spin correlations. 

In neutralino production effects from CP- violating phases can only occur if two 
different neutralinos are produced, e + e~ —> i ^ j. Each of the produced 

neutralinos has a polarization with a component normal to the production plane (SI 
EHETj. This polarization leads to asymmetries in the angular distributions of the 
decay products. 

In our process the only T^odd contribution originates from the production pro- 
cess. It is proportional to e[p e +p e -s^op^o] leading to €\p e +Pe-Pe N Pe R ] due to spin 
correlations caused by the mentioned neutralino polarization normal to the produc- 
tion plane. Applying momentum conservation pg = p^o + pi F allows one to extract 
the TN-odd triple product of the beam and the final lepton momenta [35] . 

T = (p e - x P e +) ■ P e - = (Pe-,P<+,P*-), (5) 

which projects out the CP-sensitive parts. The corresponding TN-odd asymmetry 
is then 

AT) = °(T>o)-*(T<o) (6) 

where a is the cross section for neutralino production and decay, Eqs. AS}-®- Note: 
in the case of the 3-body neutralino decay, Xi ~^ one a ls° obtains T N -odd 

contributions originating only from the decay process Therefore further TN- 

odd asymmetries can be composed that contribute also in case of same-neutralino 
pair production and can offer tools for disentangling the different phases. The CP- 
sensitive asymmetries, Eq. ([6]), depend on the charge of the leptons [13] and the 
following relations are given: 

APe-,Pe+,P e -) = -A(Pe-,P e -,P e +) 

= -A(p e -,p e -,p £ +) 

= +A(p e -,p e +,p e -), (7) 

neglecting FSI contributions. Note that a true CP-odd asymmetry, where also an 
absorptive phase from FSI or finite-width effects is automatically eliminated, can be 
defined as 

A CP = \{A-A\ (8) 

where A denotes the TN-odd asymmetry for the CP-conjugated process [68]. In our 
case this leads to a separate measurement of asymmetries for a positive {ij^j and 
negative (t^) near lepton. If A CF ^ holds than we observe a genuine CP-violating 
effect. Therefore, it is important to tag the charge of the near and far leptons in 
order to establish CP violation in the process (J2])-(|3]). 
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Figure 2: Dependence of the asymmetry A, Eq. fl6]), (a) on the phase fa (with 
</V» — 0), (b) the phase (p^ (with fa — 0), for neutralino production e + e~ — >■ X1X2 
(solid), and e + e~ — >■ X1X3 (dashed), and subsequent decay x° — > £r£n, @r Xi^f, 
(for £ = e or /i), at i/i = 500 GeV and polarized beams (P e -,P e +) = (0.8, —0.6). 
The other MSSM parameters are given in Tab. [TJ In the left panel, along the flat 
line of the asymmetry (solid) the decay X2 — * ^R^ is closed. 
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166 GeV 
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Table 1: MSSM parameters of the benchmark scenario at the electroweak scale, 
see Sec. 12.31 

2.3 Benchmark scenario 

For our full simulation study, we have chosen a benchmark scenario with the relevant 
MSSM parameters given in Tab. [TJ Since the phase of the Higgsino mass parameter 
is strongly constrained by EDM bounds, we have set it to zero. The value of the 
gaugino phase fa = 0.2tt approximately corresponds to the maximum of the CP 
asymmetries, see Fig. |2j The scenario was chosen to have an enhanced neutralino 
mixing close to a level-crossing of the neutralino states X2 an d X3 f° r fa = 0? an d of 
Xi and X2 f° r 0i = 7r ' which leads to large CP asymmetries. 

Further, we have assumed beam polarizations of (P e -, P e +) = (0.8, —0.6) enhanc- 
ing slightly the SUSY cross section and the asymmetries. At the same time, this 
choice suppresses the background from Wiy-pair production, a(e + e~ — > WW) = 
0.7 pb (compared with 7 pb for unpolarized beams), and also chargino pair produc- 
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masses m^o = 117 GeV 

Xi 

m,o = 169 GeV 

X2 

77V = 181 GeV 

X3 

m^o = 330 GeV 
m-± = 146 GeV 

Xi 

m-± = 330 GeV 

X2 





= 166 GeV 


m h 


= 280 GeV 


m fl 


= 165 GeV 




= 280 GeV 


m p 


= 268 GeV 




= 268 GeV 



cross sections 


a(e + e" -+ xlxl) 


= 244 fb 


a(e + e~ e+e^) 


= 304 fb 




a(e+e- -> xixl) 


= 243 fb 


cr(e + e- -> /ifc 


) = 97 fb 


branching ratios 


BR(x° ->> M 


= 55% 


BR(x° 2 -+ hr) 


= 45% 




BR(x§ i R £) 


= 64% 


BR(x° 3 -+ nr) 


= 36% 


asymmetries 




-9.2% 


^CHxSxg = 


7.7% 



Table 2: Masses, production cross sections, neutralino branching ratios and asymme- 
tries, Eq. ([6]), for the benchmark scenario, see Tab. (TJ calculated using the formulas 
presented in [42]. The ILC cross sections are for yfs = 500 GeV and polarized beams 
(P e -,P e +) = (0.8, —0.6). The branching ratios are summed over £ = e, fi and both 
slepton charges. 



tion a(e + e~ — > xtxj) — HO fb (410 fb, respectively). Thus, the scenario is opti- 
mized to yield large asymmetries, large cross sections and sizable neutralino branch- 
ing ratios into electrons and muons, as listed in Tab. [2j We have set A T = —250 GeV 
in the stau sector, which has low impact on the neutralino branching ratios. Also 
we have chosen the slepton masses such that x® is close in mass with the slepton t^. 
This leads to soft leptons from the decay x\ ~~ ^ £r^i as can be seen in Fig. El which 
has to be taken care of in the lepton identification described in Sec. 14.21 

In general, our analysis is relevant for scenarios with strong gaugino-higgsino 
mixing in the neutralino sector, usually leading to sizable asymmetries. In particular, 
f° r I/ 1 1 ~ \M 2 \ < 300 GeV the asymmetries can reach several percent, the neutralino 
pair-production cross sections reach more than 50 fb, and the neutralino branching 
ratios into electrons and muons are of the order of several 10%, see also Ref. [35] 
for more details. In any case, the selectron and smuon masses should fulfil m^o < 
mj> R < tu-^o, so that at least one relevant decay channel remains open. Decreasing 
the selectron mass will result in larger asymmetries and production cross sections. 
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Figure 3: Energy distributions (each normalized to 1) of the near lepton £n (solid), 
and the far lepton lp (dashed), from neutralino production e + e~ — > XiXi f° r ( a ) i = 
2, and (b) i — 3, with subsequent decay Xi ~^ ^r^n 
scenario as given in Tabs. [T] and [2J 



9 R -> X^F, 



for the benchmark 



3 Kinematic selection of signal and background 

In order to measure the CP asymmetries, Eq. ([51) , we have to separate the signal 
lepton pairs originating from X1X2 an d X1X3 production, respectively. This is essen- 
tial, since in our scenario the corresponding CP asymmetries, A(XiX2) = —9.2% and 
A(xiX,3) = 7.7%, have opposite sign. Large CP asymmetries naturally occur when 
the neutralinos are mixed states of gauginos and Higgsinos, which often implies that 
they are close in mass. In addition we need an efficient method for background 
separation. The CP-even backgrounds will reduce the asymmetries, since they con- 
tribute to the denominator, but cancel out in the numerator of the asymmetries, see 
Eq. ([6]). Among the most severe SM and SUSY background processes are W pair 
production and slepton pair production. 

Owing to a known center-of-mass energy and a well-defined initial state one may 
attempt to perform a full kinematic reconstruction of the events at the ILC. Unlike 
in the case of the LHC [24"tl2~T] , this is already possible with very short decay chains. 
Assuming that the masses of intermediate and invisible particles are known from 
other measurements, the full reconstruction can be performed even when only two 
particles are visible in the final state. 

In the following, we extend the method of Ref. [58] to reconstruct the pair of 
signal leptons from the neutralino decay, Eq. (j2J). We show that even for two neu- 
tralino states that are close in mass, here X2 an d xls? the final leptons can be correctly 
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assigned to their mother particle. A similar procedure to identify and suppress back- 
ground from W and slepton pair production is described in App.[A] see also Ref. [69|. 
Finally, we discuss how well the kinematic selection and reconstruction works at the 
MC level. 



3.1 Kinematic constraints from neutralino production 

In the center-of-mass system of neutralino pair production the momenta and energies 
are fixed 



s + m 2 Q -m 2 s + m 2 -mi A 2 (s, m 2 , m 2 ) 

771 X i X-j TP A i Xi I I M Aj /„x 
J^oO = 7= i £><>0 = 1= i PvO = 7= i 9) 

x * l^l ' 2^ ' ' v ; 

with the beam energy _E = y/s/2, the neutralino masses mg>,m^o, and A(x,y,2) = 
x 2 + y 2 + z 2 — 2(xy + xz + yz). The neutralino production is followed by the two-body 
decay chain of one of the neutralinos x° via a slepton, 

Xi -> 4 + ^v^X? + ^f + ^v, ^ = e,/i. (10) 

In our signal process, Eqs. ©-(jl]), we have x*j — Xi an d it escapes undetected. 

In the following, we assume that the near and far leptons can be distinguished 
via their different energy distributions^] For our scenario, the leptons from X1X2 
production and decay have distinct energy ranges, see Fig. |3]^a). The leptons from 
X1X3 production and decay only have a small overlap in the energy window Eg G 
[18, 38] GeV, see Fig. [31(b) . Events are discarded if both leptons happen to fall into 
this energy range. We now choose a coordinate system such that the measured 
momenta are 

P e N = \ PiN \ (0,0,1), (11) 
Pe F = \Pi F \ (sm9 NF , 0, cos9 NF ), 6 NF e[0,ir], (12) 

where O^f is the angle between the near and the far leptons. In order to fully 
reconstruct the event, the decay angles of the sleptons need to be resolved 

Pi = \Pg\ (sinfecos-B, sin&sinl?, cos6), 6g [0,7r], Be [0,27r]. (13) 

The slepton momentum, p~ = Ej — m 2 , is already fixed due to energy conservation, 
Ei = EyO — Ei N . Using also momentum conservation, p 2 = (p^ + p^) 2 , the polar 
angle, b = <X (p^, p^), can be determined 

p|q ~ P 2 e N - Pj . r _ 

coso = — — r; — j , sin = +v 1 — cos z 0. (14) 

2|P^llPll 



6 This is not needed for the determination of the asymmetry, see Eq. ijTj). but will be exploited 
for the event selection. 
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Using also momentum conservation in the slepton decay, = p^ F + p^o , a similar 
relation can be obtained for the azimuthal angle 



n 1 ( Pj + Pj ~ p|q 

COS i3 = : ; ; COS COS C/jy F i (15) 

smftsmSjvjf \ 2|p£ F ||p|-| 

and p? = E-o — m? is obtained from energy conservation E^o = E^o — Ee. T — Et-. 
The kinematic variables cos 6 and cos £> solely depend on the center-of-mass energy 
s, the lepton energies, Eg N and -E 1 ^, the angle between the leptons, p^ • p£ F , and 
finally on the contributing particle masses, m^o, m^o, and Tn^ R . Thus, there only 
remains an ambiguity for sin B, since for B £ [0, 2ir] we have sin£> = ±\/l — cos 2 B. 
This ambiguity in the azimuthal angle is irrelevant for the efficiency of the event 
selection^. 



3.2 Method of kinematic event selection 

For a given lepton pair, we apply the following kinematic selection method. The aim 
is to assign the correct origin of the lepton pair, which can be signal, e + e~ — > XiXv 
or background e + e~ — > W + W~, t^t^. Thus, we have four systems of equations, 
one for each possible production process. For each candidate event we employ the 
following kinematic selection: 

• We apply the reconstruction procedure from Sec. I3.1[ assuming X1X2 an d X1X3 
production. Thus, we calculate cos b and cos B, Eqs. (TH|) and (fT5|) . with 

m-o = m^o (m^o = m^o) for X1X2 (X1X3) production. 

• We apply the reconstruction procedure from App. |Aj assuming WW and slep- 
ton pair production. Thus, we calculate two values of y 2 , Eq. (]A.9|) . 

• The event solves the system of equations if 

I cos b\ < 1 and | cos B\ < 1 , (17) 
when neutralino production has been assumed, and 

y 2 > , (18) 

when PU/slepton production has been assumed. 
The ambiguity is related to the neutralino momentum, for which we have two possible solutions 

(|p^| sinfrcosi? \ 
±|pj| sin&sinS . (16) 
|P*J + \Pl\cosb J 

For this reason, the neutralino production plane cannot be resolved in e + e~ — > XiX^j processes 
for j = 1. For j > 2, the decay of x® could be included in order to reconstruct the production 
plane [58]. In that case larger triple product asymmetries can be studied, which include the 
neutralino momentum itself |35) . 



10 



system solved only 









y?yS 

Al A2 


y?yS 

Al A3 


^R^R 


W + W~ 


rr> 


64 k 


xlxl 


41566 


788 


64 


856 


rue procei 


74 k 


x\xi 


100 


25513 


369 


873 


200 k 


^R^R 


181 


1801 


43919 


3400 




8.8 k 


w + w~ 





13 


37 


6802 






purity 


99% 


91% 


99% 


57% 






efficiency 


65% 


34% 


22% 


77% 



Table 3: The numbers of leptonic events from the pair production of neutralinos, 
sleptons and W bosons, with their identification according to the kinematic selection 
procedure at the generator level, see Sec. 13.21 The events are simulated for our 
benchmark scenario, Tab. (TJ with an integrated luminosity of C = 500 fb" 1 and 
beam polarization (P e -,P e +) = (0.8, —0.6) for y/s = 500 GeV. 

• The event is accepted and labeled as coming from a given process only if it 
solves exactly one out of the four above mentioned systems of equations, i.e. 
it fulfills condition ffTTj) for X1X2 or X1X3 production, or condition (fTS"]) for W 
or slepton production. 

In order to demonstrate the efficiency of this procedure we perform a Monte Carlo 
simulation of X1X21 XiaIs; W + W~, and t^t^ production and their leptonic decays, 
using Whizard 1 . 96 [70]. We use the MSSM parameters, Tab. [TJ with an integrated 
luminosity of C = 500 fb -1 , and a beam polarization of (P e -, P e +) = (0.8, —0.6) with 
realistic beam spectra^. In Tab. El the results of the event selection are summarized. 
Without any additional cuts, the selection method gives an excellent separation 
between the different samples at the MC level. This method is still performing well 
after a full detector simulation, as demonstrated in the following section. 

The method can be successfully applied also for different particle mass spectra. 
In our benchmark scenario, due to the small difference between Ir and y° masses, 
a separation of the near and far leptons, and of the y° an d X3 signals was rather 
straightforward, see Fig. El However, scenarios with different m^o — mj can turn 
out to be more demanding. In order to test the applicability of the kinematic 
reconstruction, we consider a scenario with vn~ lR = 146 GeV and the other parameters 
kept as in the benchmark point. It can be regarded as the worst case scenario 
since the energies of leptons from neutralino x%i X3> an d slepton decays are all in 

8 We include ISR and beamstrahlung, which slightly degrade the number of reconstructed events 
even if the correct process is assumed for a given event. Additionally, these effects will increase 
the number of false solutions, leading to wrong assignments. 
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the 10 to 80 GeV range. Therefore, in the kinematic reconstruction, one has to 
take into account two possible assignments of the near and far leptons. This, in 
principle, could result in a reduction of the efficiency. Nevertheless, in the case 
of X1X3 production the efficiency is about 40%. It drops to 8% for X1X2 pairs, 
since they likely also solve the kinematic on-shell conditions for the X1X3 production 
process. The purity, however, remains at about 90%. In half of the cases, one can 
also correctly and unambiguously assign near and far leptons. 

4 Full detector simulation study 

The next step of our analysis is passing the generated signal events and all relevant 
SM and SUSY background events through a full ILD simulation and event recon- 
struction. After discussing the preselection cuts for the leptonic event candidates, 
we apply the kinematic selection as described in the previous section. 

4.1 Detector simulation and event reconstruction 

For the present study we have performed a full simulation of the ILD detector de- 
signed for the ILC. A detailed description of the detector concept can be found in 
Ref. [7T]. The ILD is a concept under study for a multipurpose particle detector 
with a forward-backward symmetric cylindrical geometry. It is designed for an ex- 
cellent precision in momentum and energy measurements over a large solid angle. 
The tracking system consists of a multi-layer pixel-vertex detector, surrounded by 
a system of strip and pixel detectors and a large volume time projection chamber. 
The track finding efficiency is 99.5% for momenta above 1 GeV and angles down 
to 9 = 7°, while the transverse momentum resolution is S (1/pr) ~ 2 • 10~ 5 GeV -1 . 
The SiW electromagnetic calorimeter (ECAL) is highly segmented with a trans- 
verse cell size of 5 mm x 5 mm and 20 layers. It provides an energy resolution of 
(16.6±0.1)/ a/ E(GeV) © (1.1 ±0.1)% for the measurement of electrons and photons, 
and also the steel-scintillator hadronic calorimeter is highly granular and optimized 
for Particle Flow reconstruction. The calorimeters are surrounded by a large super- 
conducting coil, creating an axial magnetic field of 3.5 Tesla. 

For the simulation of the ILD, we use the ILD_00 detector model, as implemented 
in the Geant4-based Mokka [72TT74] package. We have taken into account all active 
elements, and also cables, cooling systems, support structures and dead regions. We 
have used the radiation hard beam calorimeter (BCAL) to reject forward 77 events 
at low angles. In particular the modeling of the response of the BCAL is relevant 
for the estimation of the background from events with activity in the very forward 
regions. This background was estimated by tracking electrons to the BCAL and 
determining the probability of detection from a map of the expected energy density 
from beamstrahlung pairs [75] . 

All relevant SM background^) and SUSY signal and background events are gen- 
erated using Whizard [70], for C = 500 fb" 1 and (P e -,P e +) = (0.8,-0.6). The 

9 We consider the final states listed in Tab. [5] 
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initial selection no significant activity in BCAL 
number of all tracks iV tracks < 7 

lepton selection £ + £~ pair with £ = e, n 

| cos 6*| < 0.99, min. energy E > 3 GeV 
lower energetic £ with E < 18 GeV, or 
higher energetic £ with E > 38 GeV 
higher energetic £ with £ [15, 150] GeV 

flacop > 0-27T, fl acol > 0.27T 

final preselection p™ lss > 20 GeV 
E vis < 150 GeV 
mm < 55 GeV 



Table 4: Preselection cuts, see Sec. 14.21 for details. 



Whizard generator provides an ISR structure function that resums leading soft and 
collinear logarithms, and hard-collinear terms up to the third order [75]. The beam- 
strahlung is simulated using the Circe package [77]. After the detector simulation 
the events are reconstructed with MarlinReco [78]. We have used the Particle Flow 
concept, as it is implemented in Pandora [79] . 



4.2 Backgrounds and event preselection 

In order to clearly measure the CP-violating effects in the production of neutralinos, 
we need to have a clean sample of signal events. Otherwise the CP asymmetry would 
be reduced by the CP-even backgrounds, which enter in the denominator, see Eq. fl5J). 
We therefore apply a number of preselection cuts listed in Tab. HI to reject as much 
background as possible before applying the final selection. 



4.2.1 Initial selection 

For efficient electron and muon identification, we apply the following initial selection 
on the tracks and clusters reconstructed by the Pandora Particle Flow algorithm: 

• f 1 CAL > 0.6 and — — > 0.9 for electrons; 

Etot Ptrack 

• g° AL < 0.5 and < 0.3 (0.8) for muons, with energy E > (<) 12 GeV0 

Etot Ptrack 



ln The cut on the ratio of the total calorimeter energy and the track momentum is relaxed for 
low-energetic muons, which deposit more energy in the calorimeters. This ensures a reasonably 
high muon identification efficiency even for low-energetic muons. 
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where -Eecal is the energy measured in the electromagnetic calorimeter, E tot is 
the total measured energy in the calorimeters, and Ptrack is the measured track 
momentum in the tracking detectors. E is the energy of the Particle Flow object 
assigned by Pandora, which is derived from the track momentum in case a track is 
present or from the deposited energy in the calorimeters. We require no significant 
activity in the BCAL to reject 77 events. We select those events with less than eight 
track^ij, Arracks < 7, which efficiently removes all sorts of hadronic background. 

4.2.2 Lepton selection 

We require that two of the tracks form a pair of opposite-sign same-flavor leptons 
£ + £~, with £ = e or fi. Only electrons or muons are selected with a polar angle 
I cos^l < 0.99 and a minimum energy E > 3 GeV. There is a large contribution from 
beam induced e + e~ — > , yye + e~ — > ££e + e~ {£ = e, /x, r) background event 9^1 [7T|[75]. 
The two outgoing beam electrons are high-energetic with a small scattering angle, 
while the rest of the event forms a system of low energy and mass. If the beam 
remnants escape close to the beam pipe, and cannot be rejected by a low angle veto, 
the missing transverse momentum of the event is limited, such that the remaining 
leptons are almost back-to-back in the transverse projection (T). Thus, we apply a 
cut on their acoplanarity angle 

/ p T . p T \ 

#aco P > 0.27T with # acop = 7r - arccos I ; I T I ) ' ( 19 ) 

VlP^+l \Pe-\J 

where 8 acop = for back-to-back (180°) events. Electrons or muons from 77 induced 
rr events usually have energies below 10 GeV and can therefore be suppressed by 
exploiting that the far leptons from SUSY signal decays usually have higher energies. 
To do this and to select the signal lepton pairs £ + £~ from the neutralino X2 3 decays, 
we first use their energy distributions, see Fig. El We keep events where either 
the lower energetic lepton has E < 18 GeV, or the higher energetic lepton has 
E > 38 GeV. In addition the higher energetic lepton is required to have an energy 
EG [15, 150] GeV. Since the signal lepton pairs £ + £~ originate from the same parent 
neutralino, they follow its direction in first approximation. The lepton pairs from 
SM decay processes, and also from slepton pair decays, tend to be more back-to- 
back, since the leptons originate from different mother particles. We therefore apply 
a cut on the acollinearity angle between the leptons 

#acoi > 0.2tt with # aco i = 7r — arccos ( ^ + j— - ) . (20) 

\\Pt+ \ \Pe-\J 



11 Although we expect to have only two isolated leptons in our signal events, we do not tighten 
this cut to avoid removing signal events due to overlaid 77 events. 

12 In this study we only consider the class of 77 background events where both photons from the 
e + and e~ beam interact via a virtual fermion. Interactions where the photon fluctuates into a 
vector-meson or where the photon is highly virtual and the interaction is best described as deep in- 
elastic electron scattering on a vector-meson are not considered. Since their transverse momentum 
distributions are narrower, it is less likely that they contribute to the overall background [751 18U] . 
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Figure 4: (a) Missing transverse momentum p™ lss distribution of SM background, 
SUSY background and SUSY signal after the lepton selection, see Sec. 14.2.21 (b) 
Invariant mass mn distribution of the lepton pair after all preselection cuts except 
the cut on m«. The events are simulated for C = 500 fb" 1 , beam polarization 
(P e -,P e +) = (0.8,-0.6) at v / s = 500 GeV, and MSSM parameters for our bench- 
mark scenario, Tab. [TJ 



4.2.3 Final preselection 



In Fig. 4(a) , we show the missing transverse momentum, p™ lss , distribution of the SM 
background, the SUSY background, and the SUSY signal after the lepton selection. 
The p™ lss is calculated to balance out the sum of all reconstructed transverse particle 
momenta in an event. Our signal neutralinos \2 3? but also the background sleptons, 
decay into the lightest neutralino, which escapes detection, thus giving signatures 
with high p™ lss . However, most background lepton pairs from beam induced 77 
events have a transverse momentum typically below 10 GeV, and are removed by 
the cut pf ss > 20 GeV0 Due to the escaping neutralinos, we also expect a limited 
total visible energy E V [ S in the signal events, and we apply the cut E vis < 150 GeV. 
The visible energy is calculated as the sum of all reconstructed particle energies. 

Finally we apply a cut m« < 55 GeV on the invariant mass of the lepton pair, see 



the distribution in Fig. 4(b) , after all preselection cuts, except the cut on m«. The 
signal lepton pair from x% decays has a sharp endpoint at 51 GeV (22 GeV), 
which is also exploited for mass measurements [STHS3] . The invariant mass cut 
also removes SM backgrounds from ZZ and WW production. In Fig. 4(b) , we can 
see the invariant mass peak of one of the Z bosons decaying into two electrons or 
muons, while the other decays into a neutrino pair. The WW events contribute to 
the background if they either both decay directly into same-flavor leptons, or if one 
of them decays into a r, which in turn can complete the same-flavor lepton pair in 
the event by its subsequent decay. 

The number of remaining events after the lepton selection and the entire event 
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The spike in the 3rd bin of the p™ ISS distribution in Fig. |4(a) is due to 4 77 — > 11 events that 



have a high event weight. Due to limited CPU time and the large cross section of these events, it 
is not possible to simulate an event sample corresponding to C = 500 fb _1 . The final preselection 
cuts are chosen such that this remaining high cross section background is safely removed. 
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class 


final state 


after lepton selection 


after preselection 


signal 


X°iXi -+ Xixitt {* ^ r) 


31543 


28039 




X°iXl -> j&j&U (t ^ r) 


49084 


45966 


SUSY 


H^xIxlM (^r) 


108302 


34223 




X°lXlTT 


5147 


4076 






681 


528 


SM 


llvv 


8241 


1196 




TT 


13017 


360 




il {£ ^ r) 


24113 







qq 


1380 







77 


917355 


272 



Table 5: Number of selected events after lepton selection and final preselection, for 
C = 500 fb" 1 , (P e -,P e +) = (0.8,-0.6) at yfi = 500 GeV. The MSSM parameters 
are given in Tab. [TJ 

preselection is listed in Tab. [5J The most severe remaining SM background stems 
from WW and ZZ production, while the slepton pair production is the dominant 
SUSY background. The difference in the numbers of selected x\ and x° decays is due 
to different cross sections times branching ratios, see Tab. El and due to a reduced 
muon identification efficiency at low muon momenta, which reduces the efficiency 
for the selection of x\ decays. 

4.3 Signal identification with kinematic event selection 

In order to measure our CP asymmetry from the preselected events, we now apply 
the kinematic selection procedure, which we have described for the signal in Sec. [31 
and for the WW and 11 backgrounds in App. |Aj The kinematic selection allows us 
not only to reduce the remaining SUSY background from slepton pair production, 
but also to distinguish the lepton pairs which stem from X1X2 or X1X3 production 
and decay. This will be essential, since in our benchmark scenario, Tab. (TJ the 
corresponding CP asymmetries have roughly equal size, but opposite sign, see Fig. [2j 

For each preselected lepton pair, we require that it exclusively solves only one 
of the systems of equations, as discussed in Sec. 13.21 We reject all other events 
that solve more than one system of equations. In Tab. I4.3[ we list the number of 
preselected events that fulfill this requirement. For the lepton pairs coming from 
X1X2 decays, the final signal selection efficiency is 29%, and the total background 
contamination of the selected sample is about 8%. Lepton pairs from X1X3 decays 
reach a signal selection efficiency of 27%, while the total background contamination 
of the selected sample is about 16%. 
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class 


only x°iX°2 


only XiXl 


only £+t R 


only W+W~ 


X°iX°2 X°iX°A {* ^ r) 


18343 


615 


51 


855 


X°iX°3 -> xSxS^ (* 7^ r) 


290 


20132 


372 


635 


all SUSY background 


1153 


3055 


5626 


951 


all SM background 


87 


256 


44 


81 


purity 


92% 


84% 






efficiency 


29% 


27% 







Table 6: Number of preselected events from Tab. that fulfill the requirements of 
the kinematic selection procedure, discussed in Sec. 14.31 



4.4 Measurement of the CP asymmetries 

The CP asymmetries, Eq. can now be calculated as the difference between the 
number of events N + and N-, with the triple product T > or T < 0, respectively, 



We obtain 



with the statistical uncertainty jBT~ 



„4(p e -,p,+ ,p,-)^o = -10.2 ±1.0%, (22) 

A(p e -, Pl -,p 4 )^o = +10.7 ±1.0%, (23) 

^(p e -,p,+ ,p,-)^o = +9.3 + 1.0%, (24) 

A(p e -,p £ -,p it U^ = -8.8 + 1.0%, (25) 



<MW = (26) 

and the total number of events N = N + + AL. Exchanging the near and far leptons 
gives, within the uncertainties, the same size of the asymmetry but with opposite 
sign, see Eq. (171). However, the values of the asymmetries are different from the 
theoretical values, see Tab. (2J which is mainly due to: 

1. CP-even background events cancel in the numerator, but contribute to the 
denominator in Eq. (1211) . 

2. Events are removed by the experimental selection cuts and by the kinematic 
selection procedure, which can bias the measured asymmetry. 

CP-even backgrounds shift the asymmetry to slightly lower values, whereas the se- 
lection cuts have the opposite effect. If we assume that the background contributions 
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Figure 5: The (a 
A(p e -,p e -,p e +) m o and Afa 
The cut value used in our analysis is indicated by the dashed line. In each case 10 7 
events were generated and no detector effects are included. Statistical uncertainties 
are shown. 



are known, we obtain 



M*8 



A(p e -,P e +,P £ 
A{Pe~,P iN ,P £ +)x°x 2 

-4(p e -,P<:+,P 

•4(Pe-.P*-.P*+)x?x§ 



■11.0 ± 1.0%, 
-11.6 ±1.0%, 
-11.1 ± 1.0%, 
-10.6 ±1.0%. 



(27) 
(28) 
(29) 
(30) 



Additionally, we have studied the bias due to the selection procedure by calculating 
the asymmetry after each cut and comparing the results obtained from Whizard 
with the results obtained from Herwig++ [85H87] . Both programs show consistently 
a total shift of about 2% towards larger values after the application of the complete 
selection procedure; see also Ref. [27J. The cuts remove events with a small value 
of the asymmetry inducing an upward shift. As an example, Fig. shows the 
dependence of the asymmetry on p™ lss and # a coi, indicating also our cut value. The 
shift will be taken into account in the parameter fit described in the next section. 



5 Fit of the parameters in the neutralino sector 

In the final step of our analysis, we estimate the accuracy to determine the parame- 
ters in the neutralino sector of the MSSM. These are the six free parameters of the 
neutralino mass matrix, see App. [Bj 

\M X \, M 2 , tan/3, fa, M . (31) 

We have a number of observables at hand that can be used in the fit. These are cross 
sections, masses, and asymmetries. Masses will be measured with high precision 
using different methods [31] . For the neutralino masses we assume the uncertainties 
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as in Ref. [51], since no detailed analysis has been done for our parameter point. 
However, these uncertainties are rather conservative, since we expect that in our 
scenario a similar precision can be achieved as in [88] . In case of the cross sections, 
the uncertainty is dominated by the statistical uncertainty, 



where S and B are the signal and background contributions, respectively; see Tab.0 
Since experimentally the number of events is recorded, not the cross section itself, 
we have to take into account branching ratios for the relevant decays. These will 
depend on the masses and, in case of staus, on the stau mixing angle, cos#f. The 
stau mixing angle can be obtained from r polarization measurements in stau pair 
production [73], with an accuracy of 5% [81MJT] . 

After our procedure of the kinematic event selection, see Sec. 14.31 to disentangle 
contributions from X1X2 an d X1X3 production and decay, the background contribu- 
tions are below 15%. The small uncertainties in the beam polarizations of 0.5% |56j . 
in the luminosity, and in the SUSY masses are negligible, see also Ref. [8"TH8"3"] . For 
the CP asymmetries, we have estimated relative uncertainties of the order of 10% in 
Sec. 14.41 For the fit we take into account a bias due to cuts on the asymmetry de- 
rived from the MC simulation, as described in Sec. 14.41 Thus, the analytical value of 
the asymmetry, given in Tab. [2J is shifted accordingly. Furthermore, we use Eq. (jSJ 
to calculate the measured value of the asymmetry, which is free of FSI effects. In 
summary, we have the following set of input observables and uncertainties: 



m 9 o = 117.3 ±0.2 GeV, 
m*o = 168.5 ± 0.5 GeV, 
m^o = 180.8 ± 0.5 GeV, 

a(x\xl) x BR(x° -> Id) = 130.9 ± 1.4 fb, 

a(x\xl) x BR(x° -> Id) = 155.7 ± 1.6 fb, 

a{xlxl) x BR(x° ->• l R if = 4.8 ± 0.3 fb, 

<r(x°X°) x BR(xi -> Idf = 26.3 ± 0.7 fb, 

<7(X°X°) x BR( X ° Id) x BR(^ ->> id) = 28.9 ± 0.7 fb, 

A CP (p e -,pi N ,Pi F ) m » = +11.3% ±0.7%, 

^ CP (p e -,P^v,P£ F )x°x» = -10.9% ±0.7%. 



The uncertainties for the cross sections correspond to an integrated luminosity of 
£ = 500 fb" 1 . We perform a six dimensional x 2 fit using Minuit [92~||9"3] 



where the sum runs over the input observables Oi mentioned above, with their 
corresponding experimental uncertainties 50{. The theoretical values calculated 
using the fitted MSSM parameters, Eq. ( 131]) . are denoted by Oi. The parameter 




(33) 
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dependence of branching ratios (e.g. the stau mixing angle) is also included in the 
fit, but has negligible impact. We then obtain the following fitted values for the 
MSSM parameters: 



\ M l\ 


= 150.0 ± 0.7 GeV, 


M 2 


= 300 ± 5 GeV, 


l/'l 


= 165.0 ±0.3 GeV, 


tan j3 


= 10.0 ±1.6, 


01 


= 0.63 ±0.05, 


0,, 


= 0.0 ±0.2. 



The best estimates are obtained for the \M\\ and |/i| mass parameters, since the 
neutralino states x^ anc ^ X% are mostly composed of bino and Higgsino. The 
fourth neutralino is heavy and cannot be measured, so the limit on the wino mass 
M.2 is not as good. Also a rather large uncertainty is obtained for tan/3. However, 
if additional measurements from other sectors will be added, it should be improved 
significantly. We note that the precision obtained in this study is similar to the 
results of Ref. [S3], which uses a similar set of observablesl^l 

It is remarkable that the moduli of the phases 0i and M can also be determined 
with high precision, using the CP-even observables alone. However, only an inclu- 
sion of CP-odd asymmetries in the fit allows us to resolve the sign ambiguities of 
the phases. Without the CP-odd asymmetries in the fit we would have a twofold 
ambiguity, <fti = ±0.6, and even fourfold if M ^ 0. Thus, the triple product asym- 
metries are not only a direct test of CP violation, but are also essential to determine 
the correct values of the phases. 

6 Summary and conclusions 

We have presented the first full detector simulation study to measure SUSY CP 
phases at the ILC. We have considered CP-sensitive triple-product asymmetries in 
neutralino production e + e~ — > XiXi an d the subsequent leptonic two-body decay 
chain x® ~~ >" ^r^i ^R ~~ >* Xi^> f° r ^ = e, fi. Large asymmetries typically arise due to 
strong neutralino mixing. This causes on the one side that asymmetries for X1X2 
and X1X3 production have about the same size but opposite sign. On the other side 
the strong mixing implies two close-in-mass neutralino states, that can have a mass 
separation of the order of 10 GeV. This quasi- degeneracy would potentially pose a 
problem for the separation of both signal components. 

Therefore we have developed a kinematic selection method, to identify the lepton 
pairs from the signal events. At the Monte Carlo level, we have shown that this 

14 Thc high precision achieved in the fit calls for the inclusion of higher order corrections which 
can be in the O(20%) regime in the neutralino system, see e.g. [94j|95]. These corrections will in 
turn depend on the full parameter set of the MSSM. Therefore, the proper treatment would require 
the inclusion of observables from other sectors, in particular from the third generation of squarks, 
cf. Ref. [MlinZ]. This issue is beyond the scope of this paper, however, it should stimulate further 
studies. 
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method allows one to separate the leptons from the two signal processes X1X2 an d 
XiX®, and also to reduce the major SM and SUSY backgrounds, in particular from 
iy-pair and slepton-pair production. 

Then we have performed a detailed case study, which includes a full ILD detector 
simulation and event reconstruction. A detailed cut flow analysis has been done to 
preselect leptonic event candidates, which then have been passed to our method 
of kinematic selection. Even after the detector simulation, our method has worked 
efficiently to reduce background and separate the signal. After the full simulation 
with kinematic selection, the efficiencies of signal event selection is of the order of 
27% with a purity of about 90% of the event samples. That allows one to measure 
the asymmetry with a relative precision of about 10%. Our method of kinematic 
event reconstruction also works well in scenarios with different mass splittings of the 
neutralinos and the selectron. In the worst case scenario we found that the efficiency 
will go down to some 10%, but still with a high purity of the correctly identified 
signal sample of the order of 90%. 

We have performed a global fit of the neutralino masses, cross sections, and CP 
asymmetries to reconstruct the MSSM parameters of the neutralino sector, including 
the CP phases. The relative uncertainties of the parameters \M\\ and are below 
1%, those for M 2 about 1%, and for tan/3 and the CP phases <p±, (p^ about 10%. 
Although the moduli of the phases <fri, can also be determined by using the CP- 
even observables alone, we have shown that only an inclusion of CP-odd asymmetries 
in the fit allows us to resolve the sign ambiguities of the phases. 

To summarize, we have shown that a measurement of the neutralino sector seems 
to be feasible, including CP phases. In particular the triple product asymmetries 
are not only a direct test of CP violation, but are also essential to determine the 
correct values of the phases in the neutralino sector. 
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A Reconstruction of W and i pair production 

We consider a template process 

e + + e' ->■ A + A -> 1 + 1+ B + B, (A.l) 

where (A,B) = (£,Xi) or iW,u). In both cases B = Xi or v escapes detection. 
Since the system of the lepton pair has to obey different kinematic constraints, we 
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consider the question, whether the final lepton pair can be assigned to its mother 
production process, if the lepton momenta are measured, and the slepton and LSP 
masses are known. We follow closely Ref. [69], and define the notation 



ci = Pa- Pi 
c 2 = Pa- Pe 



1/ 2 

2 (m * 



m~ 



— ( m B — m . 



m 

4 



Ai 



(A.2) 

(A.3) 
(A.4) 
(A.5) 



h = PA- PA 

an = Pe ■ P^, a 12 = Pe ■ Pi, a 22 = Pi ■ p^- 
The momentum can be decomposed into the final lepton momenta 

Pa = t 1 p e + t 2 Pi + yp±, 
where p± = pe x pi- The expansion coefficients follow from Eqs. (IA.2j) and flA.3j) 

a 22 ci — a 12 c 2 



(A.6) 



ci 
c 2 



tlOn + t 2 d\ 2 

t\(ii 2 + t 2 a 22 



ti 
h 



Qll022 — a l2 

oiic 2 — ai 2 ci 

2 

a\\(l 22 — (1\ 2 



(A.7) 



We finally obtain, from Eqs. (1A.4j) and flA.7j) . 

^2 = (*?an + 2t 1 t 2 a u + t\a 22 ) + y 2 |pj 2 , 
b 2 - (tfau + 2t\t 2 a\_ 2 + t\a 22 ) 



y 



P-L 



(A.8) 
(A.9) 



The equation for y constitutes a condition for existence of physical solutions of the 
system, i.e. y 2 > 0, where y 2 is computed from the kinematic variables s, tjia, 



E#, Ei, and p^ ■ p^. Similar to neutralino pair production, Sec. 13. 1[ Eq. (jTojl . there 
remains a twofold ambiguity in solving the W or I system, y = ±\/y 2 . 



B Neutralino mixing 

The complex symmetric mass matrix of the neutralinos in the photino, zino, Higgsino 
basis (7, Z, H®, H®), is given by 



( M 2 s 2 w + M lC 2 w {M 2 -M 1 )s w c w 
(M 2 - Mi ) s w c w M 2 c 2 w + M lS 



\ 








m z 




m z 



\ 


-yUC 2/3 



(B.1) 



with the short hand notation for the angles sw = sin 9w, c-w = cos^iy, and s 2( g = 
sin(2/3), c 2/ 3 = cos(2/3), and the SU{2) gaugino mass parameter M 2 . The phases of 
the complex parameters Mi = |Mi|e J9i:L and fi = |/i|e^ M can lead to CP-violating 
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effects in the neutralino system. The phase of M 2 can be rotated away by a suitable 
redefinition of the fields. We diagonalize the neutralino mass matrix with a complex, 
unitary 4x4 matrix N, 

N* ■ M x o ■ N j = diag( 

with the real neutralino masses < m v o < m v o < m v o < m v o. 
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